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Abstract 

Gravitational lenses are presently playing an important role in astro- 
physics. By means of these lenses the parameters of the deflector such 
as its mass, ellipticity, etc. and Hubble's constant can be determined. 
Using C, Xforms, Mesa and Imlib a computer program to visualize 
this lens effect has been developed. This program has been applied 
to generate sequences of images of a source object and its correspond- 
ing images. It has also been used to visually test different models of 
gravitational lenses. 

1 Introduction 

Computer visualization is nowadays an important field of scientific research, 
because it often allows a better understanding of natural phenomena. 
Gravitational lenses (hereafter GL), predicted by Einstein's general relativity 
theory, deflect the light rays coming from a distant object. GL's allow the 
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formation of mutiple images of the same source object. Moreover this gravi- 
tational effect opens the possiblity of determining the mass of the deflector 
and the computation of the age of the Universe (Hubble's constant), as well 
as other lens parameters and cosmological parameters. 
An interactive computer program to visualize this gravitational effect has 
been developed. The program was written in the C programming language 
and uses the Xforms (tool kit), Mesa (graphic library) and Imlib (image 
library). The program runs on systems with Linux or Unix. An SGI version 
of this program is also available. 

In the first section the gravitational lens effect is briefly reviewed. The vi- 
sualization program is presented in the second section. In the third section 
some applications are shown: image sequences and easy visual modelling. 
Improved versions of this program can be produced and the author wishes 
to invite the interested reader to participate in this process. The author is 
preparing a website for downloading this program: 
http : //lia . ef is . ucr . ac . cr/~f rutos/ 

2 The Gravitational Lens Theory 
2.1 The Gravitational Lens Equation 

Due to the curvature of space-time, light rays coming from a distant source 
object (quasar, star or galaxy) are deflected when passing close to a lens or 
deflector (star, galaxy or cluster of galaxies). For weak gravitational fields 
the Post-Newtonian approximation is applicable. Under this approximation 
the deflection angle or Einstein's angle does not depend on the direction of 
propagation and the trayectory are approximated by straight lines. This 
gravitational light deflection is depicted in figure 1. From this figure an 
equation, the GL equation or the ray tracing equation 1 which is obeyed by 
light rays passing near a lens object can be deduced: 




(1) 



where 





2 THE GRAVITATIONAL LENS THEORY 



3 



ri D s f3 £ D s 
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and £o is the length scale on the lens surface. The angle ot(x) is called the 
Einstein's angle. The surface density of the lens is given by 

-(x) = ^ , (4) 
where S cr represents the critical density: 

c 2 D s 

A generalization of equation (0) which considers the perturbation of a ma- 
crolens (a galaxy or cluster of galaxies) is given by 

y = .M-x-a(x) , (6) 
where the matrix A4 is given by the expression 



M 



k — 7 cos (p —7 sin <p 
—7 sin 1 — k + 7 cos 1 



The components of the matrix depend upon the parameters k and 7 and 
which are respectively the dimensionless constant surface mass density, the 
dimensionless shear of the macrolens and the shear angle. The microlens (a 
star in a galaxy or a galaxy in a cluster of galaxies) is represented by means 
of the Einstein's angle a(x). 

2.2 Gravitational Lens Models 

To simulate Einstein's angle there are several models for the lens mass dis- 
tribution. The gravitational lens models are divided as follows: 

• Parametric models 

• Non-parametric models. 
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Figure 1: The Post-Newtonian approximation of the light deflection 

At the present time the program has only parametric models 1-3 . The non- 
parametric models are lately used to model gravitational lenses 4-6 . For weak 
lensing the Kaiser- Squires method can be applied to model lenses 7 ~ 10 . In 
the program these techniques have not been implemented. 
The program includes the following parametric models: 

• Chang Refsdal Model, 

• Double Plane Lens, 

• Transparent Sphere, 

• Singular Isothermal Sphere, 

• Nonsingular Isothermal Sphere, 

• Elliptical Model, 

• King Model, 

• Truncated King Model, 
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• Hubble Model, 

• De Vaucouleur Model, 

• Spiral Model, 

• Multipole Lens, 

• Rotation Lens and 

• Uniform Ring. 

The interested reader is refered to the literature for a discussion of these 
models 1-3 . 

3 The Visualization program 

3.1 Description of the program 

The program was written in C and the Mesa graphic libraries (free version 
of the Open GL) have been used. These graphic subroutines are available 
for Unix or Linux systems. The XForms Library was used to designed the 
control panel program 12 (see fig. 2). The Image Library permits to load an 
image file on the program. These libraries can be found at the following 
addresses: 

http : / /bragg . phys . uwm . edu/ xf orir 
http : / /www . mesa3d . org7 

f tp : //f tp . enlightenment . org/pub/enlightenment/ imlib/| 

The program creates a window: the Control Panel (see fig. 2). The user 

can control all item on it by clicking. A second window, the Image window, 

appears when the Image Window button is clicked on this panel (see fig. 3). 

The Help button on the panel give the user a concise program guide. 

A version of this program which employs the SGI Graphic Libraries 11 is also 

available. 

3.2 The Control Panel 

All the objects are adjusted interactively with the mouse. This control panel 
interface has the following items: 
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Figure 2: Control Panel 
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Figure 3: The Image Window 
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• Source Menu (Filled Circle, Coloured Rings and Image File) 

• Model Menu (the aforementioned models) 

• Source Positioner 

• Three counters, two for source positioning (X, Y) and one for adjusting 
the pixel resolution (N) 

• Two inputs for positioning the observed lense images (X, Y) 

• Buttons (Model Paramters, Image Window, Source, Images, Smooth, 
Grid, R Plot, R_E, File, Add, Replace, Remove, Clear, OBIP, Save, 
Help, Exit), 

• Sliders (if one chooses a model, then the corresponding sliders appear) 

• A browser for showing the observed lense images 

For a detailed discussion of these items see the manual program (to get the 
program manual send an email to the author or click the Help button on the 
control panel). 

3.3 The Image Window 

On this window the events (images, ray plot, etc.) appear. All variations of 
the parameters on the control panel are shown immediately on this window. 
To show the images of a chosen gravitational lens model the ray tracing 
method has been used. The caustic is represented by means of the ray plot 
method 1 . 

4 Applications 

In this section two of many applications of the program are discussed: 

• Image sequences 

• Easy visual modeling 

An image sequence for an elliptical model and a visual modeling of the gra- 
vitational lens 2237 + 0305 will be shown. 
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4.1 Image sequences 

An image sequence is easy to generate with this program. The structure of 
the images for different positions of the source can be investigated by means 
of such an image sequence. 

An Elliptical Lens 

A sequence of images is shown in figure 4. The sequence begins on the top 
left-hand-side panel of the sequence. The source, concentric rings, moves 
from left to right. The source is not shown, only the images. The best 
possible resolution given by the control panel was used. In this sequence the 
shear 7 and the constant surface mass density o are non zero and the shear 
angle = 0. There are two circles: the Einstein ring or lens scale and the 
core radius. The formation of a small and a large arc can be seen. In the 
middle of the sequence one can recognize the deformation of the Einstein ring 
in an ellipse. Due to the core radius an image appears in the centre, which 
fuses with an arc, in the case that the source moves away. 

4.2 Easy Visual Modeling 

Here an easy visual modeling of the gravitational lens system 2237 + 0305 
is shown. The model images are considered well-fit when they overlap most 
of the area of the observed images (represented through concentric circles), 
as visually estimated. At the present time a fitting subroutine is not imple- 
mented. When such a subroutine is included, this program would become an 
excellent tool to model gravitational lenses. 

Visual Model for the Gravitational Lens 2237 + 0305 

The system 2237 + 0305 [see fig. 5] was discovered in 1985 13 . Due to its 
unusual shape this lens is called Einstein cross. Attempts to model this lens 
have been carried out 14-20 

The data used in tables [I] and|^, was taken from Kent and Falco, and Schnei- 
der et al. 15-16 The cosmological distances of table |2| were calculated using the 
equations from Schneider et al. 15 (Hq = 100 kms -1 Mpc -1 and go — 1/2). 
Since the positions of the observed images [table |l[] are too small to be shown 
in the image window, a scale factor of 200 was chosen (the scale factor can 
be chosen by the user). Then all positions are multiplied by this scale. The 
positions are also rotated around the X axis. 
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Figure 4: An elliptical lens 
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gl2237.gif 



Figure 5: Gravitational lens 2237 + 0305 
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m2237e.gif 



Figure 6: Visual fitting of the gravitational lens 2237 + 0305 
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Table 1: Image positions for the lens 2237 + 0305 



Positions 
Object X Y 

(") (") 



2237 + 0305 A 
2237 + 0305 B 
2237 + 0305 C 
2237 + 0305 D 
Galaxy 



0.08 ±0.01 
-0.60 ±0.02 

0.70 ±0.01 
-0.77 ±0.02 

0.00 ±0.02 



-0.94 ±0.01 
0.74 ±0.01 
0.26 ±0.01 

-0.41 ±0.02 
0.00 ±0.02 



The lens galaxy is elliptical, so it is reasonable to choose an elliptical model 
in the program, although one can also try other models 15 ~ 18 . A lens model 
of N (N = 6724) point masses was used by Schenider et al. 15 . An elliptical 
King model and an elliptical de Vaucouleur model were used in Kent and 
Falco 16 . In Rix et al. a the de Vaucouleur model was used 17 . A general SIS 
model with shear was used in Wambsganss and Paczynski 18 . 
The lens major axis has a particular orientation, therefore the coordinate 
system is rotated around an angle 9, so that the elliptical potential has the 
same orientation. The angle 9 is taken as a parameter. The constant surface 
mass density and the shear are zero. The visual fitting is shown in fig. 5. In 
table [| the model parameters of the elliptical model are shown. These were 
found through variation of sliders and source positioners. 
The data for the lens scale {Einstein Radius), the core scale and the source 
position must be divided by 200. So one obtained = 0".90 , 9 C = 0".35 and 
(X, Y) = (0".14, 0".06) . The Einstein radius and core radius are determined 
by means of the equations He = oleD\ and r c = 9 c Di (7Ze = 0.48 Kpc and 
r c = 0.19 Kpc). 

From the following equation the lens mass can be found 

2 _ AGM DiDis 

The velocity dispersion is determined through an equation from Blandford 
and Kochanek 21-22 . 

The major axis angle can be calculated from 9' = 90° + 9 (9 — —24°). 
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Table 2: Parameters of the lens galaxy and the source 



Parameter 



Value 



Zl 

Distance Di 
Ellipticity e 

Angle of the mayor axis 
Angle of the bulk 
Velocity dispersion a v 

z s 

Distance D s 
Distance Di s 



0.0394 
110 Mpc 



0.57 
77° 

38° 



215km/s 

1.695 
867 Mpc 
825 Mpc 



In table |4] all values for the different models are shown: the elliptical King 
model (model A) 16 , the de Vaucouleur model (model B, the source position 
is referred to the image A) 17 , the generalized SIS model (model C) 18 and the 
elliptical model. The models A, B and C were fitted. Our elliptical model is 
not reliable, because the program does not possess a fitting subroutine, but 
one can estimate and compare with other models. 

5 Conclusions and Future Work 

From both the didactical and scientific point of view a program to visualize 
the gravitational lenses is useful. This versatile program works quickly and 
interactively with the mouse. 

With this computer program the user has a tool to visualize and to visually 
model gravitational lenses. The applications of the program, we have showed 
in this paper, are: 

• Sequences of images 

• Easy visual modeling 

The user can produce sequences of images for a chosen gravitational lens 
model. Through the variation of model parameters he or she can investigate 
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Table 3: Model parameter of the elliptical model and source position 



Parameter 




Value 


E: Lens scale 




180 


9: Rotation angle 




336° 


7: Shear 







a: Constant surface mass 


density 





0: Shear angle 




0° 


R: Source Radius 




50 


C: Core scale 




70 


e: Ellipticity 




0.54 


a: Softness 




0.05 


k: Constant Central mass 


density 


1.00 


(X, Y): Source Position 




(28.0, 12.0) 



the structure of the images. The user can also attempt to visually model ob- 
served gravitational lenses. The observed position data can be used as input 
data and the model parameters can be easily varied in order to approximate 
the observed images. So the user can quickly obtain model parameter esti- 
mations. Some observed lenses have already been modeled and the user can 
compare those results with the output from a chosen model of the control 
panel. 

Future Work 

The program can be improved by the inclusion of some additional subrou- 
tines: 

• Contour subroutine for the isochrones (time delay) 

• Light curves subroutine (dependence of brightness with time) 

• Subroutine for computing the image magnification 

• Subroutine to calculate critical curves and caustics 



Fitting subroutine 
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Model 



Parameter 


A 


B 


C 


Elliptical 


Einstein Ring 


0".90 


0".90 


0".874 


0".90 


Shear 






0.02 




Angle of the mayor axis 


66.8 


68 


66.84 


66 


Core Radius kpc 


0.10 






0.19 


Ellipticity e 


0.42 


0.3 




0.54 


Softness a 








0.05 


Central surface mass density k 








1.00 


Velocity dispersion km/s 


170 


209 




109 


Lens Mass 10 W M Q 


1.2 


1.08 


1.49 


1.15 


Source Position 










X 


-0".07 


0".159 




0".14 


Y 


-0".02 


0".877 




0".06 



• Root finder subroutine 

• Subroutine to load images of observed gravitational lenses 

• Subroutine with more complex (elliptical) models 

• Subroutine for superposition of models in different lens planes 

• Subroutine with cosmic string lens models 

• Subroutine for non-parametric reconstruction 

• Kaiser- Squires Subroutine 

The author is working on implementing some of the abovementioned im- 
provements. 
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